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We consider the generation of the scalar mode of the metric perturbations during preheating 
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stage in a two field model with the potential V(cj>, \) — 2 — I" 9 2 X • We discuss two possible 
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sources of such perturbations: a) due to the coupling between the perturbation of the matter field S\ 
and the background part of the matter field xo(t), b) due to non- linear fluctuations in a condensate 
of "particles" of the field X- Both types of the metric perturbations are assumed to be small, 
and estimated using the linear theory of the metric perturbations. We estimate analytically the 
upper limit of the amplitude of the metric perturbations for all scales in the limit of so-called broad 
resonance, and show that the large scale metric perturbations are very small, and taking them into 
|/-v ■ account does not influence the standard picture of the production of the metric perturbations in 

inflationary scenario. 
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In the modern scenario of creation of the matter in the Universe there is a stage called preheating, when the 
matter (usually imitating by a massless scalar field) is generating from vacuum fluctuations due to the effect of 
parametric resonance [1,2,3]. Namely it is assumed that the matter field x is coupled with inflaton field 0, and the 
coupling constant g is large. During inflation and right after that period the matter field has a large effective "mass" 
determined by this coupling, and each mode of the matter field with sufficiently small wavenumber evolves by a 
standard manner, oscillating with a frequency proportional to this "mass" and adiabatically decreasing amplitude. 
The inflaton starts to oscillate itself after inflation, and near the moments of time when <fi(t) « the matter field 
effectively loses its "mass" , the adiabatic approximation breaks and a possibility of resonant growth of the matter 
field amplitude and the corresponding occupation numbers of x "particles" appears. In fact, this effect may lead to 
exponential growth of the occupation numbers, and the resulting distribution of x "particles" is strongly nonthermal. 
After some moment of time £* (end of the first stage of preheating) the energy density of the x "particles" becomes 
comparable with the energy density of the inflaton, and back reaction processes become influence the dynamics of 
inflaton field </> and the Universe. It is believed that during some short period of time St after £* the x particles are 
thcrmalized by rescattering effects, and after thermalization the energy density of these "particles" evolves according 
to the standard picture of the Hot Big Bang model. Such scenario of the first stage of preheating is called the scenario 
of broad resonance [2,3]. 

As was mentioned by a number of authors (e.g [4]), the characteristic value of the growth rate of the matter field 
modes does not depend significantly on the value of wavenumber k for the modes with sufficiently small k. At first 
glance this fact gives a very interesting possibility of amplification of the matter field modes and metric perturbations 
at exponentially large scales, say, corresponding to the scale of the present horizon. In this paper we consider the 
simplest model of preheating with two scalar fields, the inflaton field <j> and the matter field x an d the potential 
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V(4>, x) = + 2 x , and show that in this model the value of metric perturbations at the scale of present horizon 
is suppressed by an extremely small factor, and therefore taking into account the stage of preheating does not lead 
to any modifications of the standard picture of generation of cosmological perturbations in inflationary scenario. We 
have two basic arguments supporting this conclusion. At first, the field x has a large effective mass before the stage 
of preheating (much larger than the Hubble parameter at the end of inflation) , and as a consequence the spectrum 



of initial field fluctuations Sxk is strongly suppressed at small k at the time ij„ of the beginning of the preheating 
stage, the rms value of the field fluctuations Sxrms oc fc 3//2 (also [5]). During preheating the field modes may grow 
exponentially fast, but the characteristic growth rate G = -p^ is constrained by some maximal value Umaxiri, where 

numerical constant [i max = — In (1 + y2) ~ 0.28 [3]. The first stage of preheating ends at the time t* ~ 504- IOOto -1 , 
and this estimate does not depend strongly on the parameters of the theory (see [3], and also the eqns. (22), (23) 
below). Thus, at the end of the first stage of preheating, the r.m.s value of the field fluctuations contains a factor 

SXrmsiQ < GSXrmsiUn) ~ e -^/2+/W,mt, = ^47^1.5(^-50)4-0.28^,-100)^ (1) 

where the N is the number of e-folds, and we choose the standard value N = 50 to represent the scale of present 
horizon. Here we express the amplitude of the field perturbation in terms of the natural Plank units. 

Secondly, the suppression factor for the metric perturbations may be even much smaller than the estimate (1) if 
one uses the standard theory. In this theory the contribution of the second field x i n t ne scalar mode of metric 
perturbations is determined by the terms containing multiplication of homogeneous "background" part of the field 
Xo(t) and the perturbed part 6x, and their derivatives (see eqns. (25,26) below). The amplitude of homogeneous 
part of the field behaves like the mode Sxk with wavenumber k = 0, and is constrained by the fact that the field xo 
cannot contribute to the total energy density during the last stage of inflation. Assuming that during last N e-folds 
the dynamics of the Universe has been controlled by the inflaton field <j>, after the end of the first stage of preheating 
the amplitude xo should contain the same factor (1). Therefore the rms amplitude of the metric perturbations 5 rms 
contains the factor (1) squared: 

A (+ \ s /O -94 /s -3(JV-50)+0.56(mt«-100) n 

u rms\ L * ) \ c e , z* 

and is suppressed by enormously small factor e~ 94 for the typical parameters of the theory. This estimate is however 
changed in a more self-consistent approach to the evolution of the perturbations during preheating. In the theory of 
preheating the role of background is effectively played by a condensate of x particles with relatively large wavenumbers 
( a typical wavenumber of such particles k cr n is always larger than the characteristic wavenumber corresponding to 
the scale of cosmological horizon during the preheating stage, see [3] and also the eq. (13)). The fluctuations of 
energy-momentum tensor of such condensate can give rise to additional metric perturbations, which are second order 
with respect to the perturbations of the matter field, and therefore cannot be obtained in the frameworks of the 
standard perturbation theory. The estimate shows that the fluctuations of the energy-momentum tensor decrease 
with scale proportional to k 3 / 2 (similar to the rms value of 5x field, see eq. (47)), and is of order of unity at the scale 
corresponding to the critical wavenumber k cr u at the time i* . Since in the large scale limit the metric fluctuations are 
of the order of the energy density fluctuations, the estimate 5 rms (t^) ~ Sxrms{t*) is a more reliable estimate in the 
more realistic approach to the dynamics of preheating, and this estimate still contains a very small number e -47 for 
the scale corresponding to the present horizon. In fact, it can be shown (see Section 3) that the metric perturbations of 
this kind take their maximal value at the scale corresponding to the horizon size at the time of the end of preheating, 
and this value is smaller than ~ I0~ 3 for the broad resonance case []. Therefore in the simplest models the transition 
from inflation to the hot Big Bang proceeds smoothly, without significant generation of the metric perturbations. 

We present our arguments in a more rigorous form in the next two Sections. In Section 2 we introduce the basic 
ideas of preheating theory, and estimate the time £* of duration of the first stage of preheating. In Section 3 we 
discuss the application of the theory of cosmological perturbations to our case of two interacting fields, and estimate 
the upper limit on the metric perturbations. 

We use below the natural system of units, and set the Plank mass M p i = \/8tt. 

II. PREHEATING IN THE REGIME OF BROAD RESONANCE 

The theory of initial stage of preheating has been developed in the paper [3], and we will closely follow this paper. 
For our purposes we need to know the evolution of both fields <f> an d X: and also the evolution of the scale factor a(t). 



1 Note, that when estimating the metric perturbations induced by the non- linear terms, we do not take into account the 
contribution of such fluctuations in the dynamics of the background model. This contribution is of order of the leading term 
at the end of preheating, and can change our estimates on a numerical factor. We believe that this factor is of order of unity, 
and cannot change our results significantly. 



As usually we divide both fields on background parts and perturbations <f> = 4>o(t) + S<fi(t,x),x — Xo(t) + Sx(t, x), 
and apply Fourier transform to the perturbed parts. We temporary neglect the influence of the metric perturbations 
in this Section, and consider this effect later on. Assuming that the contribution of the field \ m the energy density 
and the pressure is negligible, we can use the standard expression describing the evolution of the scale factor and the 
field (J)q in the theory of massive scalar field (see e.g. [6]): 
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where <j>i n = 2\ /| is a characteristic value of the field in the beginning of preheating stage, ao is an "initial" value of 

the scale factor, dimensionless time r = mt, and the preheating stage begins when approximately r = r, n w 1. 
The evolution equation for the perturbed part of the field % can be conveniently written in the form: 



X k + u? k X k = 0, 5 

where X k is the rescaled field amplitude: X k = a 3 ^ 2 5xk, ^k is the effective frequency: 

k 2 sin t 

"l = - + g 2 tin( — ) +a, 6 

a z t 

and the correction A = — (|iJ 2 + |-). Hereafter H = ^ is the expansion rate. The positive frequency solutions of 
this equation X + determine a vacuum state, and must have a form: 



(2tt) 3/2 % /2^ 



at the moments of time sufficiently close to Tj„ |j Here & — J rfiw|. The solution of the equation (5) can also be 
represented in another form by introducing two complex functions a, /?, which are constrained by the normalization 
condition \a\ — \/3\ = 1: 

(2tt) 3/2 V2^ V27J 

To reconcile this representation with (7), we should set a(r ~ r in ) = 1, and (3(t ~ T in ) = 0. Obviously, the 
representation (8) is specially convenient if the solution of eq. (5) is close to its adiabatic approximation, and 
a,/3 fa const. For the general case, the evolution of the functions a, (3 follows from the eq. (5): 

a = ±e 2l9 (3, p=± e -™a. 9 

If the adiabatic condition -^ < 1 is satisfied, the functions a, /3 are approximately constant, and no additional 
"particles" of the field \ are produced, the field oscillates with the frequency w u, and the field amplitude decays as 
a~ 3 ' 2 ~ r -1 . The adiabatic approximation breaks when the field <p is close to zero, and the time r is close to Tj = tvj 
(the integer index j must be much larger than unity for the validity of the approximate equations (3,4)). Rewriting 
the adiabatic condition near the points r = Tj, we have 

uj 2 ~ g<f> in AT 2 2q(r,) 1/2 Ar 2 ' 



2 A care should be taken when specifying the positive frequency solution for the modes with wavelengths larger than the 
horizon scale. Strictly speaking the vacuum state should be the standard Bunch-Davies vacuum state for a massive field, but 
a more accurate expression give essentially the same result. 

3 Wherever it is possible we will drop the index k in our expressions. 



where At = r — Tj, and the parameter q(r) characterizes the strength of the resonance [3]: 
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We will also use go = <z( T = -0 to parameterize our expressions. The regime of broad resonance corresponds to very 
large values of qo, see below. 

Thus when the parameter q is very large, the functions a and (3 are approximately constant between the moments 
of time Tj, but when the time is very close to tj, the functions a and j3 can be changed. The rule of change can 
be written as an iterative mapping between the functions olj-i, Pj-i and Uj, (3j corresponding to the time periods 
Tj-i < t < Tj and Tj < t < Tj+i [3]: 



a b 

b* a* ) ~ [ 8 j ~ 1 
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where a = Vl + e"™V*, and b = i e -% K2+2i8J , and 6 j = j t] dtu, tj = Tj/m, <p = argT(±*f£) + ^-{1 + 21n £), and 
k = k/k*, fc* is a characteristic value of comoving momentum: fc* = 2 1 ' 2 a(t)q 1 ' 4 m. If the wavenumber is much larger 
than a critical cutoff wavenumber: 

kent = — = \ -a(t)q 1/4 m = \ —aoqJ T l/& m, 13 

7T V 7T V 7T 

the effect of particle creation is exponentially damped. It is important to note that the bulk of the produced particles 
has values of wavenumbers close to the cutoff value (13), the corresponding wavelengths of these particles are well 
inside the cosmological horizon. The particles occupation numbers rij = \Bj\ can also be expressed by an iterative 
way: 



Uj w (1 + 2e _7rK - 2sin6t te-%K 2 y/l + e~™?)nj-i, 14 

where 9 to t = 20 J ' — (j> + arg(/3 J ) — arg(aP~ ), and the limit of large occupation numbers %_i S> 1 is assumed 
hereafter. In particular, from the eq. (14) it follows that the number of "particles" can either increase or decrease 
depending on the value of to t- The evolution of 6 tot is rather complicated in the limit of large g, but for our purposes 
it is sufficient to note that the occupation number cannot increase more than 3 + 2v2 times, and the amplitude of 
the field ~ y/rTJ cannot increase more than 1 + y2 times. It is convenient to characterize the change of the amplitude 
by the growth rate: 

1 , n j 
fij = ^— In — ^-r, 15 



and also average the growth rate over the time: 



2tt r^'- 1 ' 



^ = ^E^ 16 



From the eqns. (14 — 16) it follows that neither fj,j nor \f^ can exceed the maximal value [i max = - In (1 + \/2) ~ 0.28. 
In fact, the effective growth rate /i e ^ is about two times smaller than its maximal value. The numerical calculations 
give /j, e ff in the interval 0.1 — 0.18 with an average value of order of 0.14 for the coupling constant g in the interval 
from 0.9 -lO" 4 to 10~ 3 [3]. 

The background field xo(*) obeys the same equation (5) provided the wavenumber k is set to zero. Therefore we 
can obtain the expression for the evolution of the background field combining the positive and negative frequency 
solutions 5x+, <5x- = <5x+- Let denote the value of the field xa at the beginning of the preheating stage as Xim an d 
its time derivative as Xin ■ The values of these quantities are determined by their evolution at the previous stage of 
inflation. Assuming that the contribution of the field x i n the potential has been negligible during last N e-folds, we 
estimate the initial amplitude as Xin ~ — e~ ' 2 , and its time derivative as Xin ~ w(t ~ Ti n )xin ~ me ' . Let 
introduce a characteristic field amplitude and angle: 



L + g~ 2 \h , <^n = arctan — — . 17 



X* = \ Xin + 9 *X 
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Then 

x3/2 1/2 



Xo(t) = (2TO0) , "° X * (e^X- - e-** X+ ) 

zv2 



X^ock 3 / 2 ,^ 1 / 2 



£(— ) ( — ) (e- ie (/3*e^ -ae-^ )+e 4e (a*e^°-/3e-^ )), 18 



where the wavenumber fc is set to zero in the expressions for a and /3, and u>o = u>(t = Ti n ) ~ g. 

We also need the evolution of the field Sip, but we discuss it in the next Section. 

Since our order-of-magnitude estimates (1), (2) are exponentially sensitive to the value of the moment of time t* 
of the end of the first stage of preheating, it is very important to obtain a reliable estimate of t* . Following [3] we 
assume that the first stage of preheating ends when the vacuum expectation value for \ 2 gives a contribution to the 
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potential of order of the leading classical term m £ : 



(X 2 )(U) = -. 19 
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Let estimate the dependence of (x 2 ) on time. Using the definition of the function f}, we have 

where k cr n is given by the eq. (13). In order to calculate the integral in the eq. (20) we assume that each mode with 
a given k grows with its average growth rate /z e " w 0.14 when fc < k cr it(t), and the growth rate for the modes with 
k > kcrit{t) is zero. We have: 

[X } 12^ Me// i/2 r e 

where we approximate the value ofwaswK ^22. _ Note that our estimate is slightly different from the estimate [3] due 
to a different approximation in the calculation of the integral in (20). Substituting the eq. (21) in (19), we obtain: 

1 [2 12n^eff 1/2 r^ 

r * = mt * = 2^77 ln( V3 m2? 3/ 4 )• 22 

An approximate solution of the eq. (21) can be written as: 

r, w/i a i 4 (107 + 3.571n J F), 23 

where F — rnZ^q^ Mo 14> an d m -6 = np^i 94 — ^> A*o.i4 = *wjg ■ The time t* decreases with increase of the 
parameter q , and is bounded by some minimal and maximal values T min , T max . For example, if the mass m has 
its "canonical" value ~ 10 -6 , the parameter qg cannot exceed the maximal value q m ax ~ m~ 2 = 10 12 , and we have 
tmin ~ 55. On the other hand, the parameter <?(£*) cannot be smaller than unity for applicability of our theory. [] 
Since q ~ t -2 , we have q m in ~ 10 4 as a minimal possible value of qo. Therefore, the time of the end of the first stage 
of preheating is constrained by a condition: 55 < t* < 107 for [f^ fa 0.14 and m ~ 10~ 6 . 

III. THE METRIC PERTURBATIONS DURING PREHEATING 

Now let us take into account the metric perturbations and consider the self-consistent theory of perturbations. The 
shape of the metric perturbations depends on gauge. In many applications the gauge-independent formalism is used, 



4 Of course preheating can go on when q(r) <C 1. However, in that case the growth rate ~ q is small, and preheating is 
inefficient. In the two fields model the metric perturbations in the limit of small q have been discussed in the paper [7]. The 
similar one field model has been considered in the paper [8]. 



where all perturbed quantities are expressed in terms of their gauge- independent combinations. These combinations 
are naturally linked to some preferable coordinate system. The most commonly used preferable coordinate systems 
are so-called comoving and Newtonian coordinate systems. However, both these systems are not convenient in our 
case because the dynamical equations written in these systems have a singularity when 0o — 0. Therefore we use a 
synchronous coordinate system, and to specify the gauge we assume that our synchronous system is comoving at the 
moment of time r = Ti n . For the simplicity we also neglect the metric perturbations generated during inflation. In a 
synchronous coordinate system the metrics has the following form: 



ds 2 = dt 2 - a 2 (t)(AS% + B'^)dx a dx , 



24 



where the Greek indices run from 1 to 3, and the summation rule is used hereafter. The prime stands for differentiation 
with respect to x a . The perturbed equations of motion are (e.g. [9,10]): 



Hh = <fi 64> + xo$X 



dV , 



dV 



Sx, 
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A = -((j> 64> + xo8x), 



26 



d>ih d 2 V r 



27 



where h — 3A + AB, i,j = 1, 2, for definiteness index 1 stands for the fields cf>o, S(f>, and 2 for the fields xo, Sx- We 
neglect the terms proportional to k 2 in the eqns (25 — 27), assuming that the wavenumbers of the perturbations are 
very small. The quantities A and B are related as (e.g. [10]) 



B = f ^ J dt"(aA) 



2N 



As we discussed above the background field xo is very small, and does not contribute to the expansion law, and in 
that case the equations of motion of the fields perturbations (27) can be rewritten in a simplified manner []. The 
equation for the field Sx is reduced to the eq. (5), and the equation for the field S(f> has a form 
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J((/>,Xo,8x) = ~( 
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In this approximation the fields xoj ^X enter only in the source term J ((f), xo> $X)i anc i ^ ne formal solution of the eq. 
(29) can easily be obtained if the source term is known as a function of time J(<p(i), Xo(t), 6x(t)) = J(t)'- 
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Then one can substitute the solution (31) and the expressions for xo(t) an d Sx(t) into the eqns. (25), (26), and 
find the metric perturbations by integration. Unfortunately, this program is analytically very complicated. Even the 
numerical solution is not trivial for all possible parameters of the theory. Still, not everything is lost. At first we 
need to estimate only the upper limit on the metric perturbations, and therefore we can assume that both fields: the 
background field xo an d the perturbation Sx grow with the maximal possible rate fJ^maxm- In this approximation 



5 Obviously, when approaching the time t* the back reaction of the produced x "particles" influences the motion of our system 
much more significantly than the field xo, and comparable with the contribution of the field <j>o(t). We cannot take into 
account the back reaction without a very significant complication of our consideration, but we hope that this cannot change 
order-of-magnitude estimates. Therefore, our results should be treated as semi-qualitative only. 



at the time close to r* the growth rate is much larger than the expansion rate H = |^ <~ ^, and we can take 
into account the source term (30) only during the last Hubble epoch before the end of preheating T *~ T <C 1, t* ^> 1. 
Secondly, we should take into account only the leading terms in the expansion on powers of the small parameter 
q (t). Using the first approximation, wc neglect the terms caused by the expansion of the Universe in the equation 
of motion of the field S(f>, and write: 

, d 2 V 2 sinr 
S(p + to 5(f) = - 5\ = -2.g (pin X°X, 32 

where the term x$X can be written as: 

XS X = ** ^ 9 (-) 3/2 ^M4 + «ie" 2W + <e 2tf ), 33 
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where 



s J = (2nj + l)e**° - 2a j /3 j e- l4 ">, s J + = f3 j (a*e l<t > - /fye - '* ), si = aj-(/JV*° - Oj-e-'* ), 



and the functions <x,, /3j correspond to the period of time Tj < t < Tj+i. In the same approximation we can set 
a = a(i*) = a*. In this case during the time period Tj—i < t < tj the equation (32) is just an oscillator equation with 
a source, represented by a sum of a constant and an oscillating parts. Assuming q{t) 3> 1, it can be easily seen that 
a partial solution determined by the oscillating parts is small, and the partial solution of the eq. (32) S(f> p is close to 
a constant: 

<%, « K{-iys{, 34 

where 

_ ijg/mfg- 1 / 2 a ^/2 ^ 

V2 W (2™*) 3 / 2 ' 

and we approximate ojo ~ 5- The full solution can be found from the continuity conditions at the time moments Tj, 
and has a form: 

S(f) = K((-l) j s : ' + ^2cicosT), 36 

1=0 

where 

Now we substitute the solutions (4), (36) into the eq. (26), and by integrating the result obtain the perturbation 
of the scale factor A. Note that the term xo$X < 0(t*—4>q84>) is not taken into account. It can be said that the 
perturbation of the matter field Sx induces the metric perturbation not directly, but through the excitation of the 
inflaton perturbation S(f>. We have: 

where the function of time 

R(r) = {(-1)V sinr + \ £ q(t - ti + ^)} e - 2 <— , 38 

is constrained by the condition |i?(r)| < 0(1). The coefficient B is additionally damped with respect to A in the long 
wave limit, and is not calculated here. The vacuum expectation value for A 2 has a form: 

x 7 rf ./ 3?r 2 to t6 ./ 



where n — In (^^p) is e-folds number, corresponding to the wavenumber k, and we assume that \* ~ ~ e 3N ^ 2 . Let 
us remind that for the modes with wavelengths of order of the present horizon scale n w N « 50. The r.m.s value of 
the metric perturbation 



3-7T 2 v m ' 



A.. m ...(«) « \/^(£:) ^ me -^^ + »'+*~-, 40 



is even much smaller than the estimate (2) due to very small factor < 10 9 in the front of the exponent in the eq. 
(40). 

Thus we conclude that in our case the large-scale metric perturbation induced due to the coupling between the 
matter field perturbations 5\ and the background part of the field xo is absolutely negligible. 

There is another source of unavoidable metric perturbations induced by fluctuations of the energy density and the 
pressure of the condensate of the x "particles" . These fluctuations are non-linear, and therefore give rise to the metric 
perturbations even in the absence of the classical field xo- The calculation of this effect is even more complicated 
problem than the calculations in the linear theory, but it is possible again to roughly estimate the characteristic order 
of magnitude. Similar to the linear case, there is no hope to find some imprints of this effect at super-large scales, but 
it might play an important role right before the end of preheating, at the horizon scale A* = i? ~ 1 (£*) = |t»to _1 . The 
relatively large metric perturbations at this scale (say, with r.m.s amplitude ~ 10~ 2 -=- 10 _1 ) might lead to the copious 
production of primordial black holes, and modify the evolution of the Universe right after the end of the preheating 
stage PL An overproduction of the primordial black holes might constrain the parameters of our model. 

To characterize the energy density fluctuations at some comoving scale <r, we introduce coarse-grained energy density 
operator: 



c c.ff 



where 



e 



d 3 yG a (x-y)e(y), 41 



^(x 2 + sW) 



is the operator of the energy density of the x field, and the Gaussian window function 

G <J (f)^{2Try 3/2 a- 3 e-^. 
Consider the relative standard deviation 



The calculation of this quantity can be easily done in two limiting cases. Let introduce the wavenumber k a = a^ 1 
corresponding to the scale a, and also the wavenumber 

_ k crit _ _ [J ql /4 a m 



corresponding to the modes which have been amplified during all first stage of preheating: fco ~ k C rit{Tm), and the 
wavenumber corresponding to the horizon scale: 

u 1 \ 1 \u( \ 2a ° m 
k h {T) = a(T)H(r) = g^iyj- 



6 When calculating the abundance of primordial black holes, one should take into account that the metric perturbations 
are non-Gaussian in that case. This can lead to a modification of the standard estimates, based on assumption of Gaussian 
statistics for the perturbations. 



Let also average the quantity (42) over many oscillations of the field \. Then, if kg < k a < k cr it, the deviation S e is 
close to 1 : S e rs 1. If k a < ko , we have 

- £). 43 



e 2\/2 fco 
The eq. (43) can also be rewritten in terms of ratio k a /kh- 
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where 



^3/4 f 2 « e //) 1/4 u 1/4 

^ ~ Q 3/8 r- ~ 2 iU 3/8 r- ' 4 ° 

and T2 = jSs-. For the perturbations with wavelengths smaller than the horizon wavelength, the metric perturbation 
A is related to the energy density perturbation via Poisson equation: ^4- s=s — Se. Assuming that Se ~ S e (e), and 
also that at the end of the first stage of preheating the energy density of x "particles" influences the expansion law: 
H 2 ~ (e) , we have: 

A(k a )~(^p±) 2 S e = D(^p±) 1/2 , 46 

for kh(r*) < k a < fco- For super- horizon perturbations the dependence of the perturbations on k should be the same 
as a dependence of a source of the perturbations, and 

A{k a ) ~ S e = D(-^-) . 47 

Thus, the characteristic metric perturbation induced by the fluctuations of the energy density takes its maximal value 
at the horizon scale, and this value is smaller than ~ 10~ 3 if qo > 10 4 . 

Similar to the case of linear perturbation, we did not take into account the contribution of \ 'particles' in the 
dynamics of background model when calculating the amplitude of the metric perturbations induced by the non-linear 
terms. However, this contribution can change the parameters of the background model (such as, e.g. the expansion 
rate) only on factor of order two, and we believe this modification cannot change significantly our estimates f\. 
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